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Abstract 


This paper will discuss the application of system identification techniques and robust control strategies to a proton exchange membrane fuel-cell 
system. The fuel-cell system’s dynamic behaviour is influenced by many factors, such as the reaction mechanism, pressure, flow-rate, composition 
and temperature change, and is inherently non-linear and time varying. From a system point of view, however, the system can be modelled as a 
two-input, two-output linear time-invariant system whose inputs are hydrogen and air flow rates, and whose outputs are cell voltage and current. 
On the other hand, the system’s non-linearities and time-varying characteristics can be regarded as system uncertainties and disturbances that are 
treated by the designed robust controllers. This paper is comprised of three parts. First, system identification techniques were adopted to model 
the system’s transfer functions. Second, the H,, robust control strategies were applied to stabilise the system. Finally, the system’s stability and 
performance were compromised by introducing weighting functions to the controller’s design. From the experimental results, the designed Hx 


robust controllers were deemed effective. 
© 2006 Elsevier B.V. All rights reserved. 
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1. Introduction 


As the levels of available fossil fuel decrease, alternative 
energy resources gain more and more attention. Among them, 
the fuel cell is an important candidate to replace traditional fuel 
because of its favourable features, such as low operational tem- 
perature, fast power response, high power density, low noise 
pollution, high system efficiency and environmental friendli- 
ness. Until recently, fuel cells have been applied to many systems 
such as bikes, motorcycles, vehicles, boats, etc. [14,15]. As a 
result, the control of fuel cells to increase system efficiency and 
stability is becoming more and more important. The traditional 
methodology for controlling fuel-cell systems is mainly linear 
control algorithms such as PI based diagonal controllers. How- 
ever, these methods cannot fully satisfy the complex load change 
required for many applications, and may cause instability prob- 
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lems. Therefore, more elegant control strategies, such as robust 
control, should be applied. 

Rodatz et al. [1] presented a dynamic model of air supply 
and developed an LQG regulator for a proton exchange mem- 
brane (PEM) fuel-cell system. From those results, performance 
improvements with the PI controller were noted. Moreover, the 
pressure trace was successfully decoupled from the mass flow 
trace, while a faster response time was also achieved. Pukrush- 
pan et al. [2] developed an observer-based feedback controller to 
protect the fuel-cell stack from oxygen starvation during changes 
in current commands, while the linear quadratic technique was 
employed based on the linearised state-space model. Sedghisi- 
garchi and Feliachi [3] developed an Hæ controller to regulate 
the system’s output voltage under small load variations. In their 
simulations, the output voltage variation was kept below 5% by 
controlling the hydrogen flow rate. Caux et al. [4] were able 
to control the air supply of a proton exchange membrane fuel- 
cell (PEMFC) system under varying current loadings. Those 
simulations proposed a species balance model to maintain con- 
stant pressure on the cathode (oxygen) compartment and to 
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follow a desired air flow-rate. Takeuchi et al. [5] emphasized 
the power management of fuel-cell systems interconnected with 
a utility grid, while successfully introducing typical topolo- 
gies of fuel-cell based distributed generation systems. Jurado 
and Saenz [6] developed an adaptive controller of a fuel-cell 
micro-turbine hybrid power plant. Considering the models’ 
time-varying dynamics and the plants’ disturbances, the adap- 
tive controller designed at a fixed operating point was utilised 
to stabilise the system under different operating conditions. 

The fuel-cell system demonstrated in this paper is designed 
and manufactured by Chung Shan Institute of Science and 
Technology (CSIST) and assembled by DELTA Electronics. It 
consists of 15 cells that are connected with pre-treated Nafion® 
112 membranes under a hot press process for optimum condi- 
tions. The platinum loading is about 0.2 mg cm7? at the anode 
and 0.4 mg cm~? at the cathode. The cells are electrically con- 
nected in series. The system’s inputs are hydrogen and air, while 
its outputs are cell voltage and current. 

The fuel-cell system’s dynamics is non-linear and time vary- 
ing. When the current load is changed, the voltage varies from 8.5 
to 12 V as shown in Fig. 1, using the on-board controller. Because 
a steady power source is important for electrical equipment, 
however, robust control methodologies are selected in order to 
provide a steady voltage or current when the operating condi- 
tions change. Robust control is well known for its capability of 
dealing with system uncertainties and disturbances [8—10,12]. 
Furthermore, the implementation of robust controllers is rather 
simple because the parameters are constant. Compared to other 
advanced control strategies such as adaptive or gain-scheduling 
control [13], in which the parameters need to be updated accord- 
ing to the system’s responses, robust controllers are much easier 
to implement. 

This paper is arranged in the following sections. In Section 
2, the fuel-cell system’s dynamics is introduced, and further 
simplified as an SISO linear system by the system identifica- 
tion techniques, in which the un-modelled dynamics is treated 
as system variations. Section 3 discusses robust control strate- 
gies and their design to achieve the systems’ maximum stability 
bounds. Section 4 introduces weighting functions and explains 


Voltage variation as load current changed: 0A 732A 373A 374A > 5A 
T T T 


4 T T T T T T 


voltage (V) 


9Ļ 


8.5 1 1 f 1 1 f 1 1 1 
0 100 200 300 400 500 600 700 800 900 1000 


time (sec) 


Fig. 1. The voltage variation as the load current is changed. 


their improvement of closed-loop systems’ overall performance. 
Finally, Section 5 discusses some conclusions that have been 
drawn from this analysis. 


2. System identification of the fuel-cell system 


In this section, the fuel-cell system will be introduced and 
simplified as a linear model. The system transfer functions are 
obtained by a dynamic signal analyser, HP 35670A, and will 
then be used for the controller design in the following sections. 


2.1. System description 


The fuel-cell stack is designed and manufactured by CSIST 
and assembled by DELTA Electronics. It consists of 15 cells 
with the active area of 50cm? on each on them. The cells are 
connected with a pre-treated membrane — Nafion® 112 — by 
hot press for optimum conditions. Platinum loading is about 
0.2mgcm~? at anode and 0.4 mg cm~? at cathode. The cells 
are electrically connected in series, with rated and peak power 
outputs of 117 W at 9 V and 124 W at 7.8 V, respectively. The 
maximum efficiency of the fuel cell stack is 37% (LHV) under 
dry H2/air and humidification-free conditions. 

The dynamics of the fuel-cell system is influenced by many 
factors, including the diffusion dynamic, the Nernst equation, 
proton concentration dynamics and cathode kinetics as illus- 
trated in Fig. 2: 


Diffusion dynamic : Rohm = Rref + ot (T — Tef), (1) 


Nernst equation : 


dE? RT 1/2 
E = Esef + qr f — Tref) + koe In( Py, Po, ), (2) 
Proton concentration dynamics : 
OC y+ OCy+ Cyt 1 + yt ’P 
u = : (3) 
ot ot TH+ TH+ 


Pio [H+ jr 
Cathode kinetics : 1 = b In { a P (: = +) } l 
(4) 


From this point of view, the system can be represented as 
the MIMO block diagram shown in Fig. 3 with the following 
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Fig. 2. The dynamics of the fuel cell system. 


706 F-C. Wang et al. / Journal of Power Sources 164 (2007) 704-712 


Nair 
4] e 4 


Nin 


Fig. 3. The block diagrams of the fuel cell system. 


relations [11]: 
jee = G1 Nair + G3 Np3, (5) 
Veen = G2Nair + G4Np, + Ricell, (6) 


where G1—Gy, represent the input—output relations of the system. 
Since the dynamics of the fuel-cell system is non-linear and 
time varying, the system of (5) and (6) is inherently non-linear 
and time varying as well. By fixing the hydrogen input and 
current output, the system is further simplified as a SISO system 
in the following: 


Veet = G2Nair- (7) 
2.2. System identification 


The non-linear system G2 of (7) can be approximated by 
linear models at various operating points, given the input and 
output signals of the system. When considering a linear system 
with the sinusoidal input r=a cos(wt), as shown in Fig. 4, the 
output y are measured at N sampled points such that the system 
transfer functions can be calculated by the system identification 
method described as follows [7]: 


N 
I= TDO cos(wt), (8) 
t=1 
1 N 
Is = 5) VO sinod, (9) 
t=1 
a VIN) + P(N 
eyes ee (10) 
a/2 
LN) 


én Go) = argG y(e) = — arctan 


I(N)’ an 
where |Ĝ y(e!”)| and @y (jw) are the gain and phase of the system 
at frequency øw. 

To model the system transfer function more efficiently, a 
Dynamic Signal Analyzer (DSA) HP 35670A is utilized. At first, 
a swept-sine analogue voltage signal, which is generated from 
the analyser, is sent to both of the air pump and channel-one of the 


rG 


Fig. 4. A linear system. 


Table 1 

The operating conditions of the system 

Load (A) Input 

H2 Air 

2 1.2(LPM), 7 psi 50 mHz-100 Hz sine wave 
(level: 2 V, offset: 4 V) 

3 1.2(LPM), 7 psi 50 mHz-100 Hz sine wave 
(level: 2 V, offset: 4 V) 

4 1.2(LPM), 7 psi 50 mHz-—100 Hz sine wave 


(level: 2 V, offset: 5 V) 


DSA to obtain the mathematical models of the system. Since the 
current of the input signal is in the scale of 1 mA, while the driv- 
ing current of the air-pump measures approximately 500 mA, a 
power amplifier circuit is constructed to actuate the air pump. 
Under different operating conditions, such as those illustrated 
in Table 1, the output voltage of the fuel-cell system can be 
connected to channel-two of the DSA for system identification. 
Furthermore, to capture the system dynamics, the frequency 
range of input signals is set at between 0.05 and 100 Hz with 
400 sampled frequencies. At each sampled frequency, the sys- 
tem gain and phase are calculated from (8) to (11), allowing the 
system transfer functions to be estimated from the frequency 
response diagrams. The experimental set-up of the system is 
illustrated in Fig. 5. 

The frequency responses diagrams (Bode plots) of the fuel- 
cell system, under different loading conditions (2 A, 3 A and 
4 A), are illustrated in Fig. 6. Here, the approximate linear trans- 
fer functions of first-, second- and third-order are selected to 
fit the frequency responses appropriately. It is noted that the 
air-pump dynamics was also modelled in the transfer func- 
tions during the identification procedures. It is arranged in this 
way because the air-pump dynamics cannot be singled out in 
the control process. Lower-order models of the system can be 
used if the high frequency responses are insignificant for the 
systems. 

Three experiments were executed for each operation condi- 
tions to consider the system variations: the resulting transfer 
functions are illustrated in Tables 2-5. Those linearised models 
will be used for a discussion of robust controller design in the 
following sections. 
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Fig. 5. The experimental setup for system identification. 
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Fig. 6. The Bode plots of the system under three different operation conditions. 


Table 2 
Identification of the first-order transfer functions 
No. Load (A) 
2 3 4 
i a 0.0010s + 0.2630 ee 0.0170s + 0.6192 bes 0.00145 + 0.3254 
eS 540.2741 2 s+ 0.6501 3 = s+ 0.5244 
0.00235 + 0.2502 0.01815 + 0.6427 0.00235 + 0.3895 
2 2 = = Aaa Gy = —— G32 = ———___ 
s + 0.3521 s + 0.5913 s + 0.8171 
5 ae 0.0050s + 0.3294 Bree 0.01298 + 0.5959 ba 0.00015 + 0.3386 
1B s+ 0.4557 3 5 +.0.5501 3 5 +.0.7354 
Table 3 
Identification of the second-order transfer functions 
No. Load (A) 
2 3 4 
j P 0.00085? + 0.2372s + 0.0198 t 0.00865? + 0.4856s + 0.5118 ő 0.00405? + 0.3774s + 0.9233 
1 = 324 0.31875 + 0.0438 a s2 + 1.2193s + 0.6338 31 52 4 3.45658 + 1.4539 
> ” 0.00295? + 0.3073s + 0.1879 0.00525? + 0.4289s + 2.6402 - 0.00435? + 0.4298s + 0.5780 
12 2 F 1.23245 + 0.2176 = 52+ 4.3915 + 2.5485 2 42.411 + 1.1781 
3 y 0.0039s? + 0.3160s + 0.3955 ` 0.00105? + 0.2181s + 9.5332 7 0.00515? + 0.4499s + 0.9693 
B= S24 1.61825 + 0.5636 23 2 + 16.1856s + 8.9642 3 52 + 3.90355 + 1.9979 
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Table 4 
Identification of the third-order transfer functions 
No. Load (A) 
2 3 4A 


0.0026s? + 0.28065? + 0.59435 + 0.1943 


0.0049s? + 0.41635? + 2.64875 + 0.4103 


0.00235? + 0.32635? + 5.47545 + 11.1666 


1 Gy = Gn = = 
ss S? + 2.84725? + 1.30945 + 0.2734 s s? + 5.378182 + 5.17965 + 2.0519 L O S E 15.22019? + 43.56565 + 17.4807 
2 G 0.00365? + 0.31295? + 0.8772s + 0.9834 0.001553 + 0.23308? + 6.46095 + 2.5733 0.00315? + 0.41625? + 2.9779s + 2.7421 
12 B + 3.13725? + 5.23405 + 1.2132 = s3 + 10.9154s? + 9.97085 + 2.6322 z S + 8.113152 + 13.38345 + 5.6468 
3 Ga = 0.0017s? + 0.25565? + 2.9615s + 0.4103 Gn = 0.00025? + 0.04925? + 12.7968s + 0.2916 _ 0.00035? + 0.12195? + 10.6176s + 33.8400 
a S? + 9.733852 + 5.36685 + 0.6996 A s? + 21.9285s? + 12.1025s + 0.1491 3 3 22.87095? + 126.56675 + 70.0607 
Table 5 
The gap metric and stability margin of the first-order systems (with minimization of the maximum gap and the corresponding stability bound in boxes) 
Gi Gi2 G13 Gazı G2 G23 G31 G32 G33 b(G, K) 
Gi 0.0000 0.1465 0.1391 0.2829 0.3076 0.2837 0.2081 0.3150 0.3274 0.9280 
Gi2 0.1465 0.0000 0.0805 0.2887 0.3202 0.2992 0.0789 0.1727 0.1862 0.9530 
G43 0.1391 0.0805 0.0000 0.2115 0.2451 0.2257 0.0707 0.1808 0.1932 0.9522 
Gr 0.2829 0.2887 0.2115 0.0000 0.0664 0.0643 0.2288 0.3117 0.3238 0.9316 
Gn 0.3076 0.3202 0.2451 0.0664 0.0000 0.0295 0.2738 0.3730 0.3859 0.9193 
G23 0.2837 0.2992 0.2257 0.0643 0.0295 0.0000 0.2671 0.3716 0.3846 0.9186 
G31 0.2081 0.0789 0.0707 0.2288 0.2738 0.2671 0.0000 0.1110 0.1241 0.9619 
G32 0.3150 0.1727 0.1808 0.3117 0.3730 0.3716 0.1110 0.0000 0.0314 0.9756 
G33 0.3274 0.1862 0.1932 0.3238 0.3859 0.3846 0.1241 0.0314 0.0000 0.9768 
Max. 0.3274 0.3202 (0.2451 0.3238 0.3859 0.3846 0.2738 0.3730 0.3859 
Table 6 
The gap metric and stability margin of the second-order systems (with minimization of the maximum gap and the corresponding stability bound in boxes) 
Gi Gi2 G13 Gz G22 G23 G31 G32 G33 b(G, K) 
Gi 0.0000 0.2844 0.1869 0.3292 0.3696 0.3826 0.1471 0.2407 0.2452 0.9406 
Gi2 0.2844 0.0000 0.1082 0.3248 0.3323 0.3096 0.1463 0.2537 0.2582 0.9459 
G43 0.1869 0.1082 0.0000 0.2465 0.2488 0.2218 0.0469 0.1555 0.1600 0.9527 
Gr 0.3292 0.3248 0.2465 0.0000 0.1240 0.1373 0.2649 0.2792 0.3020 0.9283 
G22 0.3696 0.3323 0.2488 0.1240 0.0000 0.0372 0.2731 0.3401 0.3445 0.9126 
G23 0.3826 0.3096 0.2218 0.1373 0.0372 0.0000 0.2569 0.3524 0.3568 0.9135 
G31 0.1471 0.1463 0.0469 0.2649 0.2731 0.2569 0.0000 0.1095 0.1141 0.9622 
G32 0.2407 0.2537 0.1555 0.2792 0.3401 0.3524 0.1095 0.0000 0.0320 0.9760 
G33 0.2452 0.2582 0.1600 0.3020 0.3445 0.3568 0.1141 0.0320 0.0000 0.9776 
Max. 0.3826 0.3323 (0.2488 0.3292 0.3696 0.3826 0.2731 0.3524 0.3568 
Table 7 
The gap metric and stability margin of the third-order systems (with minimization of the maximum gap and the corresponding stability bound in boxes) 
Gi G12 G13 Gazı G2 G23 G31 G32 G33 b(G,K) 
Gi 0.0000 0.0907 0.1062 0.3274 0.3266 0.4627 0.1281 0.2219 0.2164 0.9412 
Gi2 0.0907 0.0000 0.1509 0.3238 0.3351 0.4059 0.1134 0.2277 0.2299 0.9488 
G43 0.1062 0.1509 0.0000 0.2453 0.2528 0.5382 0.0503 0.1407 0.1436 0.9549 
Gr 0.3274 0.3238 0.2453 0.0000 0.1074 0.4183 0.2624 0.2784 0.2998 0.9273 
G22 0.3266 0.3351 0.2528 0.1074 0.0000 0.3189 0.2796 0.3289 0.3403 913] 
G23 0.4627 0.4059 0.5382 0.4183 0.3189 0.0000 0.5054 0.6021 0.6039 0.9052 
G31 0.1281 0.1134 0.0503 0.2624 0.2796 0.5054 0.0000 0.1168 0.1188 0.9618 
G32 0.2219 0.2277 0.1407 0.2784 0.3289 0.6021 0.1168 0.0000 0.0359 0.9758 
G33 0.2164 0.2299 0.1436 0.2998 0.3403 0.6039 0.1188 0.0359 0.0000 0.9769 
Max. 0.4627 0.4059 0.5382 0.4183 (0.3403 0.6039 0.5054 0.6021 0.6039 


E-C. Wang et al. / Journal of Power Sources 164 (2007) 704-712 709 


3. Robust stability design 


In this section, robust control strategies will be introduced. 
From the analyses of gap metrics and coprime factorisation, 
robust controllers are designed to provide the maximum stability 
bound for the fuel-cell system. The resulting controllers are then 
verified from both simulations and experiments. 


3.1. Principles of robust control 


Robust controllers can be designed for the identified linear 
models by using robust control strategies, such as Small Gain 
Theorem, perturbations of coprime factors, stability margin and 
gap metric [8—10]. A system is said to be well-posed if all closed- 
loop transfer matrices are well-defined and proper. Furthermore, 
a closed-loop system is said to be internally stable if all the 
system states are asymptotically stable, i.e., the states go to zero 
from all initial state conditions when the input signals are zero. 


Theorem 1 (Small Gain Theorem [8]). Suppose M € RH% and 
let y>0. Then the interconnected system shown in Fig. 7 is 
well-posed and internally stable for all A(s) € RH% with 


(a) ||Alloo < Wy if and only if ||M(5)||o0 < y; 
(b) ||Alloo < 1/y tfand only if ||M(s)||o0 < y, 
where ||G||oo is the œ norm of system G. 


Suppose that a nominal plant Go can be expressed as Go = MN}, 
where (1) M,N € RH, and (2) MM“ + NN“ = Iœ. This is called a 
normalised left coprime factorisation of Go. Now, let us suppose 
that a perturbed system G4 with the block diagram of Fig. 8 is 
expressed as: 


Ga =(M+ Ay) l(N + Ay), (12) 


Fig. 8. The block diagram of G4 =(M + Am)! (N+ An) with feedback con- 
troller K. 


with ||[Ayw, An]lloo<ée, Am, An E RHoo, the system transfer 
functions can be simplified as follows: 


el-f 


K -1 

I (-— GK) [IG]o. (13) 
Hence, from Theorem 1, the closed-loop system remains inter- 
nally stable for all ||[Ajy, An ]|loo <€ if and only if 


K 
| | d (CEH a 


(1— GKY!M o 


= 


m | = 


’ 


[0,6] 
, K L1 
or equivalently I (I— GK) [IG] 


[0,6] 


(14) 


Defined the stability margin b(G, K) of the closed-loop system 
as follows: 


—1 
l (5) 


(oe) 


K -1 
b(G, K)= | l | (— GK) [IG] 


from Theorem 1 the closed-loop system is internally stable for 
all ||[Am, An]| loo <€ if and only if b(G, K) > e. 

It is further noted that the coprime factorisation of a system 
is not unique. That is, there is more than one expression for Go 
and Ga. Therefore, the gap between two systems Go and GA 
can be defined as: 


Definition 2 (Gap metric [8]). The smallest value of ||[Ay, 
An]lloo which perturbs Gg into G4, is called the gap between 
Go and G4, and is denoted as ô(Go, Ga). 


Thus, b(Go, K) gives the radius (in terms of the distance in the 
gap metric) of the largest ball of plants stabilised by controller 
K. Therefore, the design goal is to derive a suitable controller 
K from a nominal plant Go, such that all perturbed plants G; 
located inside the gap 6(Go, Gj) < € will satisfy b(Go, K) > £ and 
the closed-loop system remains internally stable. 


3.2. Robust controller design 


To perform robust controller design, several nominal plants 
are selected as follows: 
0.0050s + 0.3294 


First-ord tem, G =, 16 
irst-order system, Go_1(s) s £ 0.4557 (16) 


Second-order system : 


_ 0.0039s* + 0.031605 + 0.3955 


G = , 17 
0-2(s) s2 + 1.61825 + 0.5636 o 
Third-order system : 
0.00155? + 0.233052 + 6.4609s + 2.5733 
Go_3(s) = 3 2 
s3 + 10.9154s? + 9.97088 + 2.6322 
(18) 
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Fig. 9. The experimental setup of the PEMFC control system. 


DAQ and 


Simulink Platform 


NI SCB68 Connector 


The choices of the nominal plants Go(s) are based on the calcula- 
tion of system gaps between the nominal plants and the perturbed 
plants such that the maximum gap is minimized as: 


min max 6(Go, Gi). 
Go 


For instance, for the first-order system, G13 is chosen as the 
nominal plant because the maximum gap between G13 and 
other plants is 0.2451—the minimum for all systems (shown 
in Table 5). 

By applying robust control to the nominal plants, the resulting 
H controllers are: 
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—0.448352 — 4.65025 — 1.5751 


s2 + 11.4022s + 4.2596 
(21) 


Third-order system : K3(s) = 


The stability bounds of the controllers are b(G, K)=0.9522, 
0.9527 and 0.9131, respectively (Tables 5—7). It must be noted 
that the maximum gaps between the nominal plants and the per- 
turbed plants are always less than the stability bound, such that 
internal stability can be achieved for all system variations by 
using the designed controllers. 


3.3. Experimental results 


The designed controllers were implemented to verify the 
effects. For the experiments, a closed-loop control system of 
the PEM fuel-cell system was developed by using the National 
Instrument’s (NI) Data Acquisition (DAQ) card and Simulink. 
This was accomplished by transferring the control signal through 
the output port of the DAQ card and by measuring the volt- 
age of the PEMFC system to the input port of the DAQ 
card. A step-down circuit was constructed to divide the mea- 
sured fuel-cell voltage by 2 to protect the DAQ card since 
the allowable voltage range of the DAQ card is from —10 
to +10 V, while the maximum voltage of the PEM fuel-cell 
system is about 14 V. The experimental set-up of the PEM 
fuel-cell control system is illustrated in Fig. 9. The output 
responses of the closed-loop systems are shown in Fig. 10, 


First-order system, K1(s) = —0.3263, (19) where the controllers (19)-(21) were implemented to con- 
trol the fuel-cell system under different loading conditions 
—0.3234s — 0.3776 i i 
Second-order system : K2(s) = , (20) when the reference voltage was set at 9.5 M First of all, it 
s + 1.2079 is noted that the closed-loop systems remain stable under 
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Fig. 10. The output voltage responses of (a) K; (s) for 2 A loads, (b) K2(s) for 3 A loads, (c) K3(s) for 4 A loads, and (d) K1 (s), (e) K2(s), (£) K3(s) for current setting 


3A>2A—>4A. 
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Fig. 11. The robust controller design with weighting functions. 


system variations. Secondly, the controlled voltage of the air- 
pump was 5-6 V. Finally, the output voltage was approximately 
10-11 V, while the reference (or desired) voltage was set at 
9.5 V. 

From the experimental results, it is shown that the closed- 
loop systems are stabilised by the designed robust controllers. 
However, it is also noted that the system performance is not 
satisfactory because the controllers were designed to achieve 
the maximum stability bound, such that system performance 
was relatively ignored. For example, in the experiments the ref- 
erence voltage was set at 9.5 V, but the output voltage varied 
from 10 to 11 V. Therefore, it is necessary to consider the over- 
all performance criteria by adding weighting functions into the 
controller design. These issues are discussed in the following 
section. 
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4. Robust performance design 


The controller design discussed in Section 3 focused on the 
stability of the closed-loop systems: the design criterion was 
to discover the maximum stability bound to cope with system 
uncertainties. However, the system performance was relatively 
ignored, as shown in the system responses of Fig. 10. Therefore, 
in this section, weighting functions will be introduced to improve 
the system performance. It will be shown from the results that 
the system performance is significantly improved by sacrificing 
a little stability bound. 


4.1. Weighting function selection 


There were steady-state errors between the reference and the 
system outputs, as illustrated in Fig. 10. To track the reference 
command, an integral was added to the weighting functions. 
First, a weighting function W(s) = (s + 1)/s was added to the con- 
troller design, as illustrated in Fig. 11. For the third-order system 
with the nominal plant Go_3(s) of (18), the robust controller was 
found as: 


1.17315? + 12.91205? + 12.6523s + 3.2806 
s3 + 12.012452 + 14.1096 + 3.8368  ” 


with a maximum stability bound b(Go 3, K4) =0.6499. 
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(22) 
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Fig. 12. The output voltage responses of the fuel-cell system with K4(s) under loadings: (a) 2 A, (b) 3 A, (c) 4A and (d) 2 A > 3 A > 4A. 
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4.2. Experimental results 


The experimental results using the designed robust controller 
K4(s) are shown in Fig. 12, where the reference voltage is 
set at 9.5V while the loadings are set at 2A, 3A, 4A and 
2A—>3A-— 4A, respectively. It is shown that the closed-loop 
system is internally stable for system perturbations, and that the 
steady-state error is zero because of the integral factor of the 
weighting function. It is noted the use of the weighting func- 
tion can be regarded as a compromise between the stability 
and performance. As the system performance is improved using 
the controller K4(s), the relative stability bound of the closed- 
loop systems is decreased from b(Go.3, K3)=0.9131 to b(Go.3, 
K4)=0.6499. That is, the system robustness is reduced. 


5. Conclusion 


This paper has proposed system identification techniques and 
robust control methods for a PEM fuel-cell system. From system 
identification, the fuel-cell system was modelled as an MIMO 
system, which was further simplified to an SISO system by fix- 
ing the hydrogen input and the current output. By selecting the 
nominal plants, the system variations were regarded as system 
uncertainties and disturbances that were treated by the proposed 
robust controllers. At first, a robust controller was designed to 
achieve the closed-loop systems’ maximum stability margins. 
Then, weighting functions were introduced to improve the sys- 
tem’s overall performance. The experiments have shown that 
the designed robust controllers can achieve good system perfor- 
mance and stability. Although this paper has been limited to the 
SISO robust controller design only, it is important to note that 


the principles of developing an MIMO robust controller will be 
similar and could be applied widely. 
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